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Abstract 

We present a systematic method for proving nonterminating basic 
hypergeometric identities. Assume that k is the summation index. By 
setting a parameter x to xq n , we may find a recurrence relation of the 
summation by using the g-Zeilberger algorithm. This method applies 
to almost all nonterminating basic hypergeometric summation formu- 
las in the book of Gasper and Rahman. Furthermore, by compar- 
ing the recursions and the limit values, we may verify many classical 
transformation formulas, including the Sears-Carlitz transformation, 
transformations of the very-well-poised 8^7 series, the Rogers-Fine 
identity, and the limiting case of Watson's formula that implies the 
Rogers- Ramanuj an identities. 
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1 Introduction and Notations 



This paper is devoted to develop a systematic method for proving nontermi- 
nating basic hypergeometric series summation and transformation formulas. 
The idea of finding recurrence relations and proving basic hypergeometric 
series identities by iteration has been used very often, see [10-12,25]. How- 
ever, there does not seem to exist a systematic method within the scope 
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of computer algebra for proving nonterminating hypergeometric summation 
and transformation formulas. One obstacle lies in the infinity of the summa- 
tion ranges. In this paper, we find that the g-Zeilberger algorithm can be 
used as a mechanism for proving many basic hypergeometric summation and 
transformation formulas. 

Wilf and Zeilberger developed an algorithmic proof theory for identities 
on hypergeometric series and basic hypergeometric series [42,47,48]. For 
the purpose of this paper, we are concerned with the g-Zeilberger algorithm. 
Koornwinder [38], Paule and Riese [41], and Boing and Koepf [18] further 
studied algorithmic proofs of basic hypergeometric identities. Most of the 
theory and implementations are restricted to the case of terminating iden- 
tities. For nonterminating identities, Gessel [28] and Koornwinder [39] pro- 
vided computer proofs of Gauss' summation formula and Saalschutz' sum- 
mation formula by means of a combination of Zeilberger's algorithm and 
asymptotic estimates. 

Our approach can be described as follows. Let 

oo 

f(a,...,c) = ^2t k (a, ...,c) 

k=0 

be a hypergeometric series with parameters a, . . . , c. We aim to find a recur- 
rence relation of the form 

p Q (a, . . . , c)f(a, . . . , c) + px(a, . . . , c)f(aq, . . . , cq) H 

+ p d (a,...,c)f(aq d ,...,cq d ) = 0, (1.1) 

where d is a positive integer and Po, ■ ■ ■ ,Pd are polynomials. To this end, we 
try to find polynomials p , . . . ,p d and a sequence (g ,gi,...) such that 

p (a, . . .,c)t k (a, . . . ,c) +pi(a, . . .,c)t k (aq, . . . , cq) H 

+ p d (a, c)t k (aq d , . . . , cq d ) = g k+1 - g k . (1.2) 

Assume that go = lim g k = 0. Then (jl.lj) follows immediately by summing 

k^oo 

over k in (|1.2}1 . 

The main idea of this paper is to use the g-Zeilberger algorithm [18, 
38,41,47] to solve Pi and g k . The bridge to the g-Zeilberger algorithm is 
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the introduction of a new variable n by setting the parameters a, . . . , c to 
aq 11 , . . . , cq n . Then the summand tk{aq n , ■ ■ ■ , cq n ) becomes a bivariate q- 
hypergeometric term. Applying the g-Zeilberger algorithm, we can always 
obtain an equation of form (jl.2j) . 



When the recurrence relation (jl.lj) is of first order (i.e., d = 1) or involves 
only two terms, f(a, . . . , c) equals its limit value lim f(aq N , . . . , cq N ) mul- 

jV— »oo 

tiplied by an infinite product. Using this approach, we can derive almost all 
nonterminating summation formulas listed in the appendix of the book of 
Gasper and Rahman [27], including bilateral series formulas. 

When the recurrence relation involves at least three terms, we show that 
f(a,...,c) is determined uniquely by the recurrence relation and its limit 
value under suitable convergence conditions (Theorem 13. lj) . Therefore, to 
prove an identity, it suffices to verify that both sides satisfy the same recur- 
rence relation and that the identity holds for the limiting case. Using this 
method, we can prove many classical transformation formulas. 

Let us introduce some basic notation. The set of integers and nonnegative 
integers are denoted by Z and N, respectively. Throughout the paper, q is a 
fixed nonzero complex number with \q\ < 1. 

The g-shifted factorial is defined for any complex parameter a by 

(a; q)oo 



(a; q)oo = Y[(i 

For brevity, we write 

Oi, 



aq ), and (a; q) n 



q)n = (ar, q), 



(aq n ; q) c 



Vn e z. 



where n is an integer or infinity. Furthermore, the basic hypergeometric 
series are defined by 



ax, . . . , a 7 
h,...,b s 



q,z 



E 

k=0 



(a 1 ,...,a r ;q) k z k 
(bi,...,b s ;q) k (q;q) k 



-i)V 



s— r+1 



and the bilateral basic hypergeometric series are defined by 



r?Ps 



a i 



6i, . . .,b s 



q,z 



y> (ai, . . . , a r ; g)k ^ 
u-^ fa, ■ ■ ■ A; q)k 



■l) k q^ 
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2 Summation Formulas 



In this section, we present a method of proving nonterminating basic hyper- 
geometric identities by using the g-Zeilberger algorithm. Given a term of the 
form 

(a' l ,...,a' u ,q) 00 (a 1 ,...,a r ]q) k d Q) k 

(b' 1 ,...,b' v ;q) 00 (&!, . . . , b s ; q) k 
by setting some parameters a, . . . , c to aq n , . . . , cq n , we get a bivariate q- 
hypergeometric term t k {aq n , . . . ,cq n ) on n and k. By the g-Zeilberger al- 
gorithm, we obtain a bivariate g-hypergeometric term g n>k and polynomials 
Pi{q n , a, . . . , c) which are independent of k such that 

Po(q n , a, • • • , c)t k (aq n 1 . . . , cq n ) + p x (q n , a, . . . , c)t fc (ag n+1 , . . . , cq n+1 ) + ■■■ 
+ p d (q n , a, . . . , c)t k (aq n+d , . . . , cg n+d ) = g nMl - # n , fc . (2.1) 

Suppose that g ,o = lim g ,fc = 0. By setting n = in (|2.1|) and sum- 
ming over k, we derive a recurrence relation of form (jl.ljl for f(a,...,c) = 

oo 

]T t fc (a, . . . ,c). 

fc=0 

When the recursion (jl.ljl involves only two terms, say f(a,...,c) and 
f(aq d , . . . , cg d ), we have 

f(a, . . . , c) = lim /(ag dAr , . . . , cO • lim ff ( - Pd{ ^r ' ' . 

Therefore, the evaluation of f(a, . . . , c) becomes the evaluation of its limit 
value lim f(aq dN , . . . , cq dN ), which is much simpler and is usually an infinite 

N—*oo 

product. 

2.1 Unilateral Summations 

We now present an example to show how to obtain an infinite product ex- 
pression from an infinite summation. 

Example. The g-binomial theorem: 



The g-binomial theorem was derived by Cauchy [20], Jacobi [35] and 
Heine [31]. Heine's proof consists of using series manipulations to derive the 
recurrence relation 

(l-z)f(a,z) = (l-az)f(a,qz). (2.2) 

Gasper [26] provided another proof using a recurrence relation respect to the 
parameter a: 

f(a, z) = (1- az)f(aq, z). 

Our computer generated proof is similar to Heine's proof. The recurrence 
relation generated by the g-Zeilberger algorithm turns out to be (J2.2j) . Let 
Uk(z) be the summand 

/ \ 0;<?)fc k 
u k (z) = - — — z 

and u n) k = Uk(zq n ). By the g-Zeilberger algorithm, we obtain 

(-azq n + l)u n+ i,k + {zq n - l)«n,fc = 9n,k+i ~ 9n,k, (2.3) 
where g n ^ = (1 — q k )u ny k- Denote the left hand side of ()2.3|) by t n ^. Then 

oo 

}t n k = -g n>0 + lim g ntk = 0, Vn > 0, 

k=0 

implying that 

f(a, zq n ) = lf^.f[a, zq n+1 ), Vn > 0. 
1 — zq n 



Hence 



„. , l — az., . 1 - az 1 - azq 2 

f(a,z) = f(a,zq) = - f{a,zq) 

1 — z 1 — z \ — zq 

(az;q) N N (az^q)^ 

hm — ^ — f(a,zq ^ - 



N->ce (z;q) N ' (z;q)oo ' 

where lim f(a,zq N ) = 1 holds by Tannery's Theorem (see [45, p. 292] or 
[17]): 
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Tannery's Theorem. Suppose s(n) = J2k>o fk( n ) 1S a convergent series for 
each n. If there exists a convergent series ^ k>0 Mk such that \fk(n)\ < 
then 

oo 

lim s(n) =} lim fk(n). 

k=0 

The following summation formulas (most of which come from the ap- 
pendix of [27]) can be obtained by the above method. We only list the 
recursions and the limit values ^lim f(aq dN , . . . , cq dN ). 
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Name 


Summation 


Recurrence Relation 


Limit 
value 




00 k 

k=0 


f(z) — — f(za) 
J \ z ) — i- z J K^H) 


1 


/7-PYnnn put" 1 ^ 1 




f( z ) — (1 1 z \ f( za ) 


1 


Lebesgue 
[8, p. 21] 


(x;q) kq y > 
(<?;9)k 

fc=0 


f( r \ — n — rn) f(ra 2 ) 

J \ x ) — y 1 • LL i)j y x q ) 


(-a- a) 


Generalization 
of Lebesgue 
[6,22] 


oo ( k + 1 ) 
(a,b;q) k q\ 2 > 

(q;q)k(abq;q 2 ) k 


f(a,b) = 

(l-a q )(l-bq) n, b 2\ 
(l-abq)(l-abq 3 )-> v"^ ' U( * > 


(-91 ?)oo 


101 


, \ a c] 

i0i „ 5 5, - 
I c a\ 


f( c ) = l - c/a f(cq) 


1 


q- Gauss 


201 


a, 6 c 
c 




m = { i^:^f{cq) 


1 


g-Kummer 
(Bailey- 
Daum) sum 


201 


a, 6 g 
_ ag/6 




/(«) = 

(l-aq 2 /b 2 )(l-aq) ,/ 2 \ 
(1-aq 2 /b)(l-aq/b) J ^ > 


(-959)00 
(-q/b;q)oo 


A g-analogue 
of Bailey's 

2^l(-l) 


202 


a ' q/a h -, q ,-b 
-q,b 






1 


A g-analogue 
of Gauss' 

2^l(-l) 


202 


a 2 ,b 2 

1 1 " 

abq 2 , — abq 2 " 


f(a, b) = 

{1 -a 2 q)(l-b 2 q) f( 
(l-a 2 b 2 q)(l-a 2 b 2 q3)J ^ U( U 


9)oo 


g-Dixon sum 


/ a 2 ,—qa,b,c 
4<P3 [ -a,a 2 q/b,a 2 q/c ! 

9. a 


/(«) = 

(l-a 2 q/bc)(l-a 2 q 2 /bc)(l+aq) 
{l-a 2 q/b)(l-a 2 q 2 /b)(l-aq/bc) 
(l-aq/b)(l-aq/c)(l-a 2 q) ft \ 
(l-a 2 q/c)(l-a 2 q 2 /c) ■> ^1 


1 
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Here are more examples. 
A g-analogue of Watson's ?> F 2 sum: 

/i, qpL^, —q^, a 2 , b 2 q, c, —c, —abq/c 



f(a,c) = 807 



/12 , — /J2 , —bcq/a, —ac/b, —abq, abq, c 2 ' ^' ab 



where \i = —abc. 

By the g-Kummer sum, we have 



lim f(aq ,cq ) = 2 <f>i 

N— >oo 



b 2 q,— abq/c c 



—bcq/a 

By computation, one derives that 



ab 



(~g; q)oo(b 2 q 2 ,c 2 q/a 2 ;q 2 ) 
(c/ab, -bcq/a; q)^ 



/(«,«) = j 1 + ± g/p + «y a - y - c/t) cg) 



(1 - c 2 g)(l - abq)(l + a6g)(l + acq/b)(l + ac/b)' 



Thus, we have 
f(a,c) 



—abcq, —c/b, c/b, —q; q)oo(a 2 q, b 2 q 2 , c 2 q/a 2 ; q 2 
(abq, —abq, —ac/b, c/ab, —bcq/a; q) OD (c 2 q; q 2 ) c 



A g-analogue of Whipple's 3 F 2 sum: 

-c, g(-c)2, -g(-c)3, a, g/a, c, -d, -g/d _ 



/(c) = 807 



^— c) 2 , — (— c) 2 , —cq/a, —ac — q, cq/d, cd ' 



q,c 



By computation we have the recurrence relation 
(1 + cg 2 )(l + c)(l - cg 2 /W)(l - acg/d) 



/(c) 



(1 + acq){\ - cdq){l + cq 2 /a){l - cq 2 /d) 



(1 - acd)(l - cdq/a)(l + eg) 2 2 
X (1 + ac)(l - cd)(l + cg/a)(l - cq/dy^ ] ' 



Since /(0) = 1, we obtain 



The sum of a very-well-poised 6 5 series: 

f(a) = 6 05 



a,qa2 , —qaz ,b,c,d _ aq 
a2 , — aa ,aq/b,aq/c,aq/d ' ' 6crf 



By the g-Zeilberger algorithm, we find 

(1 — aq/cd)(l — aq/bc)(l — aq/bd)(l — aq) 



/(«) 



'1 — aq/bcd)(l — aq/b){l — aq/c){l — aq/d) 
Since /(0) = 1, we obtain 

(ag, aq/bc, aq/bd, aq/cd; q)^ 



f(aq). 



M = 



(aq/b, aq/c, aq/d, aq/bcd; q) c 



2.2 Two- Term Summation Formulas 

Many classical two-term nonterminating summation formulas can be dealt 
with by using the same method as single summation formulas. It turns out 
that for many two-term summation formulas, the two summand share the 
same recurrence relation. Moreover, the boundary values lim g ,fc for the two 

summands cancel out. So we still obtain homogeneous recurrence relations 
which lead to infinite products. We give three examples from the appendix 
of [27], and present a detailed proof for the first example. 

1. A nonterminating form of the g-Vandermonde sum: 

f(a,b,c) = 



2<Pl 



a, b 



;q,q 



+ 



(q/c,a,b; q) c 



-2<Pi 



(c/q,aq/c, bq/c; q) 
Since lim f(aq N ,bq N ,cq N ) does not exist, we consider 

N^oo 

g(a, b, c) = /(a, 6, c) / (q/ c; q)^. 

Let 



aq/c, bq/c ' 
q 2 /c ' g ' g 



u 



u 



(1) 

n,k 

(2) 
n,k 



(aq n ,bq n ;q) k k 
0. i 



(q/cq n ;q)oc (cq n ,q;q) k 



(aq n , bq n ; q) c 



{aq/c, bq/c; q) k ^ k 



(cq n /q, aq/c, bq/c; q)^ (q 2 /cq n , q; q) k 
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We have 



(abq^ - cR^ fc + cuZ = - g", i = 1, 2 



,(') 



where 



o (1) 
o (2) 



u 



n.ki 



c(l-abq 2n+k )(l-q k )^ {l) 
q k (l - aq n ){l - bq n ) 
c(cq n — q k+1 )(l — q k ) 



q k+1 (l - aq n )(l - bq n ) 



u 



(2) 
n.k' 



Noting that = = and 

,(!) 



lim 9nl = ~ , lim & 



(af+\bf+\qU 
fc-^oo-"-- fc-*oo""»~ V(V c 9 n > c n+1 >0;0)oo' 

we get g(a, 6, c) = (1 — abq/c)g(aq, bq, cq). Since 

lim ^(ag^, 6^, cg^) = + - — —— — = - — —— - 

N-*oo (aq/c,bq/c; q)oo {aq/c,bq/c;q) c 



we get 



f(a,b,c) = (q/c; q) 00 g{a 1 b,c) 



(q/c, abq/c; q) c 



(aq/c, bq/c; q) c 

2. A nonterminating form of the g-Saalschiiz sum: 



/(c) = 302 



a, b, c 



Q,Q 



+ 



(q/e,a,b, c, qf/e-q)^ 
'e/q^aq/eM/e.cq/eJ^q)^ 

aq/e, bq/e, cq/e 



3V2 



<l 2 /e,qf/e 



where / = abcq/e. 

By computation, we have 

(1 — bcq/e)(l — acq/e) 



(1 — cq/e) (l — abcq/e) 



f(cq), 



and by (I2~3J) 



lim f(cq 

N—>oo 



(q/e, abq/e- q) c 
(aq/e^q/e)^ 
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Thus we get 



f(c) 



(bcq/e, acq/e, q/e, abq/e; q) c 



(cq/ e, abcq/ e, aqj e, bqj e; q)^ ' 
3. Bailey's nonterminating extension of Jackson's 8 </> 7 sum: 



(2.5) 



f(a, b) = 8 07 



a, qaz , — qaz , b, c, d, e, / 



a*, -as, aq/b, aq/c, aq/d, aq/e,aq/f 

__ b (aq, c, d, e, f, bqj a, bqj c, bqj d, bqj e, bqj f; g)^ 

a (aq/b, aq/c, aq/d, aq/e, aq/ f, bc/a, bd/a, be/a, bf /a, b 2 q/a; q)oo 

. b 2 /a, qba~^ , —qba~^ ,b, bc/a, bd/a, be/a, bf ja 

ba~^ , —ba~^ ,bq/a,bq/c,bq/d,bq/e,bq/ f ' ' ^ 

where / = a 2 q/bcde. 

By computation, we have 

^1 — aq)(l — aq/cd)(l — aq/ce)(l — aq/de 



f(a,b) 



[1 — aq/cde)(l — aq/c)(l — aq/d)(l — aq/e) 



f(aq,bq), 



and by (j2~Kjl . 



,. ft n u tvn (b/a,aq/cf,aq/df,aq/ef;q) c 
hm f(aq ,bq ) = — — — — — 

n->oc [aq/ J, bc/a, bd/a, be/a; q)^ 



Finally, we have 



f(a,b) 



(aq, aq/cd, aq/ce, aq/de, b/a, aq/cf, aq/df, aq/ef; q) c 
(aq/cde, aq/c, aq/d, aq/e, aq/ f, bc/a, bd/a, be/a; q)^ 



2.3 Bilateral Summations 

Bilateral summations ( [27, Chapter 5]) can also be dealt with by using the 
g-Zeilberger algorithm approach. We need the following special requirement 
for the recurrence relation (j2.1j) : 



lim g n k = lim g n>h = 0. 

k— >— 00 k— >oo 
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Here are several examples. 
1. Jacobi's triple product: 



fc=— oo 



This well-known identity is due to Jacobi [34] (see [5, p. 12]). Cauchy [20] 
gave a simple proof using the (/-binomial theorem. For other proofs, see 
Andrews [1], Ewell [24], Joichi and Stanton [36]. 

We give a g-Zeilberger style proof of its semi-finite form [21]: 



f(m) = £ 



OO 11 !■» 



k=—oo 



[q m+1 \q)k 



m > 0. 



Let u m ,k be the summand. Applying the g-Zeilberger algorithm, we obtain 

zu m +i,k — (q m+1 + z)(l — q m+1 )u mtk = g m ,k+i — 9m,k, 
where g mM = (1 - q m+1 )q m+1 u m ^ k . Since 

lim g m k = lim g m , k = 0, 

k—>co k — ► — oo 



/(m + !) = (! + q m+1 /z)(l - q m+1 )f(m). 



we have 
It follows that 

oo oo 

Eq{%) z k = lim u m>k = lim f(m) 

= /(0)(g, -q/z;q)oc = (~z,q,-q/z; q)^. 



k=—oo 



k=—oo 



2. Ramanujan's sum: 



fib) = 



(q.b/a.az.q/az-.q)^ 
(b,q/a,z,b/az;q)oo ' 



|z|, |6/az| < 1. 



This formula is due to Ramanujan. Andrews [3,4], Hahn [29], Jackson [33], 
Ismail [32], Andrews and Askey [9], and Berndt [16] have found different 
proofs. 
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The proof of Andrews and Askey [9] is based on the following recursion: 

l-b/a 



fib) 



-f(bq). 



(2.6) 



(l-b)(l-b/az)' 

Instead of using series manipulations, we derive the recursion ()2.6|) by using 
the g-Zeilberger algorithm. Let u n)k = T&rh z k . Then 



Q>q n ;q) k ' 



z(bq n - a)u n+ i >k + (az - bq n )(l - bq n )u n)k = g n ,k+i ~ 9n,k, 



(2.7) 



where 



g n ,k = (1 - bq n )bq n ■ u n , k . 

Notice that (|2.7|) holds for any k G Z. Furthermore, when \z\ < 1 and 
\b/az\ < 1, we have 

lim g U)k = (1 - bq n )bq n ■ lim u n)k = 0. 

k— >±oo k— >dzoo 

Summing over 6 Z on both sides of ()2.7j1 . we immediately get 1)2.6)1 . im- 
plying that 



(b, b/az;q) c 



By the g-binomial theorem 

oo 

/(?) = E 

Therefore, 



fc=0 



{b.b/az^q)^ [q/a^q)^ 
3. A well-poised 2 ^2 series: 

/(M= 2 ^r 6 ' c 

By computation, we have 

(1 - ag/6c)(l - ag 2 /6c) 



ag 

ag/6, ag/c ' <? ' 6c 



(az;g)oo 

(^;g)oo ' 

_ (b/a,q,q/az,az; g)^ 
(b, b/az, q/a, z; q)^ ' 

\aq/bc\ < 1. 



(1 + aq/bc)(l + aq 2 /bc) 

{l-aq 2 /b 2 ){l-aq 2 /c 2 ) 



[l-q/b)(l-q/c)(l-aq/b)(l-aq/c) 



f(b/q,c/q). 
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By Jacobi's triple product identity, we obtain 

oo 

lim f(b/q N ,c/q N ) = £ q k \-a) k = (q 2 , qa, q/a; gV 
Thus, we get 



fc=— oo 



f(b,c) 



(aq/bc; q)oo(aq 2 /b 2 , aq 2 / c 2 , q 2 , qa, q/ a; q 2 ) c 



{-aq/bc, q/b, q/c, aq/b, aq/c; q) 
4. Bailey's sum of a well-poised 3 ^ 3 : 
f{b, c, d) = 3^3 



b, c, d q 



q/b, q/c, q/d ' Q ' bed _ 
We notice that applying the g-Zeilberger algorithm directly to 

(b, c,d; q) k 



(q/b, q/c, q/d; q) h 



bed) ' 



does not give a simple relation. Using an idea of Paule [40] of symmetrizing 
a bilateral summation, we replace A; by —A; to get a summation 



3^3 



b, c, d 



Q 

q/b, q/c, q/d ' ' bed 



Now we apply the g-Zeilberger algorithm to the average of the above sum- 
mands: 

l + q k (b,c,d;q) k / _q_ \ k 
2 (q/b,q/c,q/d;q) k \bcdJ 

and obtain that 

(1 - q/bc){\ - q 2 /bc){\ - q/bd){\ - q 2 /bd) 



f{b,c,d) = 



x 



(l-q/b)(l-q/c)(l-q/d) 

(l-q/cd)(l-q 2 /cd) 



(1 - q/bcd)(l - q 2 /bcd)(l - q 3 /bcd) 
By Jacobi's triple product identity, we have 



f(b/q,c/q,d/q). 



lim f(b/q N ,c/q N ,d/q N )= £ q^ (-q) k = (q\ q, q 2 ; q'U = (q; q) c 



k=— oo 
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So we get 

f{b, c, d) 



(q/b, q/c,q/d,q/bcd; q)^ 
5. A basic bilateral analogue of Dixon's sum: 



f{b,c,d) = 4^4 



— qa, b,c,d qa 



_ -a,a 2 q/b,a 2 q/c,a 2 q/d ' <? ' bed 
By computation, we get 

(1 - a 2 q/bc){l - a 2 q 2 /bc){l - a 2 q/bd)(l - a 2 q 2 /bd) 



f(b,c,d) 



(1 - a 3 q/bcd)(l - a 3 q 2 /bcd){\ - a 3 q 3 /bcd) 
(l-a 2 q/cd)(l-a 2 q 2 /cd) 
X (l-q/b)(l-q/c)(l-q/d) 



Hence, 



b c = (a 2 q/bc, a 2 q/bd, a 2 q/cd, aq/b, aq/c, aq/d; q)^ _ ^ g 
(a 3 q/bcd, q/b, q/c, q/d, a 2 q/b, a 2 q/c, a 2 q/d; q)^ 

where 

5(a) = lim f{b/q»,c/q»,d/q»)= £ t^g 3 © (-ga 3 ) fc - 

Especially, replacing 6, c, d in (|2.8j) by —a,c/q N ,d/q N and taking the limit 

N — > oo, we get 

lim /(-a, c/<A d/q N ) = . S (a). 

n^oo {-q/a, -aq; q)oo 

By Jacobi's triple product identity, we have 

oo 

lim f(-a,c/q N ,d/q N ) = V q"\-a 2 ) k = (q 2 , a 2 q, q/a 2 ; q 2 )^, 

N— >oo z — ' 

fc=— oo 

which implies that 



(q,a 2 q,q/a 2 ;q) c 

b(a) — 



(aq,q/a; q) c 
15 



Therefore, we obtain 

(a 2 q/bc, a 2 q/bd, a 2 q/cd, aq/b, aq/c, aq/d; q) c 



f(b,c,d) = 



(a 3 q/bcd, q/b, q/c, q/d, a 2 q/b, a 2 q/c, a 2 q/d; q)^ 

{q,a 2 q,q/a 2 ;q) c 



x 



(aq,q/a; q)^ 

6. Bailey's very-well-poised 6 ^ 6 series: 

, , qah,-qah,b,c,d,e qa 2 

f(b,c,d,e) =6^6 .1 i ../.,_../_ ;?>■ 



a2 , — aa , aq/b, aq/c, aq/d, aq/e ' ' bcde 

This identity is due to Bailey [15]. Other proofs have been given by Slater 
and Lakin [44], Andrews [7], Chen and Liu [23], Schlosser [43], and Jouhet 
and Schlosser [37]. Askey and Ismail [14] gave a simple proof using 6 5 sum 
and an argument based on analytic continuation. Askey [13] also showed 
that it can be obtained from a simple difference equation and Ramanujan's 
\ipi sum. 

Using our computational approach, we obtain 

(1 - aq/bc)(l - aq 2 /bc){l - aq/bd){l - aq 2 /bd) 



fib, c, d,e) 



x 



x 



(1 — aq/b)(l — aq/c)(l — aq/d)(l — aq/e) 
(1 — aq/be)(l — aq 2 /be)(l — aq/cd){\ — aq 2 /cd) 

(l-q/b)(l-q/c)(l-q/d)(l-q/e) 
(1 — aq/ce)(l — aq 2 /ce)(l — aq/de)(l — aq 2 /de) 



[1 — a 2 q/bcde)(l — a 2 q 2 /bcde) (I — a 2 q 3 /bcde)(l — a 2 q A /bcde) 

x f(b/q, c/q,d/q,e/q). 

By Jacobi's triple product identity, we have 



lim f(b/q N ,c/q N ,d/q N ,e/q N ) = - J- £ (1 - a<t*)q<$(qa 2 ) k 

fc=— oo 

= - V q^X-af = (q,aq,q/a;q) c 

1 — a L — ' 

k=— oo 

Hence we get 

(aq/bc, aq/bd, aq/be, aq/cd, aq/ce, aq/de, q, aq, q/a; q)^ 



fib, c,d,e) 



(aq/b, aq/c, aq/d, aq/e, q/b, q/c, q/d, q/e, qa 2 /bcde; q) c 
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3 Transformation Formulas 

In this section, we show that many classical transformation formulas of 
nonterminating basic hypergeometric series can be proved by using the q- 
Zeilberger algorithm. The basic idea is to find the same recurrence relation 
and limit value of two summations f(a, . . . , c) and g(a, . . . , c). Suppose we 
have obtained a recurrence relation of second order or higher order of the form 
(II. ip for both f(a, . . . , c) and g(a, . . . ,c), then the following theorem ensures 
that f(a, . . . , c) and g(a, . . . , c) must be equal as long as lim f(aq N , . . . , cq N ) 

coincides with the limit lim g(aq N , . . . , cq ). 

N^oo 

Theorem 3.1 Let f(z) be a continuous function defined on the disc \z\ < r 
and d > 2 be an integer. Suppose that we have a recurrence relation 

f(z) = ai (z)f{zq) + a 2 (z)f(zq 2 ) + ■■■ + a d (z)f(zq d ). (3.1) 

For i = l,...,d, we denote aj(0) by Wi. Suppose that there exists a real 
number M > such that 

\ai(z) — Wi\ < M\z\, 1 < % < d, 

and 

\w d \ + \wd-i + Wd\ H h \w 2 H h Wd\ < 1, 

Wi + w 2 H \-w d = l- 

Then f(z) is uniquely determined by /(0) and the functions di(z). 
Proof. By the recurrence relation (|3.1j) . we have 

d 

f(z) = J2^f(zq^), 

8=1 

where = Oi(z) and 

f 4+i = ^(zg^ 1 )^ + 4 +1) , 1 < z < d, 
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(3.2) 



Let 



X(x) = - £ £ 



.r 



i=2 



3=i 



By the assumption, A(l) > 0. Hence we may choose a real number p such 
that \q\ < p < 1 and A(p) > 0, namely, 



E^ 



j=2 



Let 



A = max{|4 1) |,...,|4 1) |}, 

A' = max {A ■ dMr/X(p), \A? - A®\/p, \A^ - 4\\/p d 
B = dMr/p d ~ 2 + A'p/A. 

We will use induction on n to show that 

\A^\<A-(-B;p) n , 
\Ag ] -4-il<^' -P n -(-B;p) n . 



(3.3) 
(3.4) 



By definition, the inequalities (|3.3|) and ()3.4|) hold for n = 1, . . . , d. Suppose 
n > d and the inequalities hold for 1,2, ... ,n. From ()3.2|) it follows that 



\A (1) I 

l^n+ll 



i=l 

d 



i=2 \ j=j / 
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By the inductive hypotheses, it follows that 
d 



< ^Mrlf+^l-A-i-Bip^-i + A-i-Bw)* 



i=l 



n+2-i 



i=2 



E 



2=1 



W, 



< A(l + dMr/p d ~ 2 ■ p n )(-B; p) n + A ■ (-B; p) n £ 



P 



n+2-i 



i=2 



E 



w , 



< A(l + (dMr/p d - 2 + A'p/A) ■ p n ) (-B; p), 
= A-{-B;p) n+1 . 

Similarly, by the inductive assumptions we have 

d 



i=l 



i=2 \ j=i J 



< A ■ dMr/p d - l -p n +\-B-p) n + A ■ (-B;p) n J2 



P 



n+2-i 



i=2 



E 



j=i 



W , 



= A [AdMr / (A p d ~ l ) + 1 - \{p)/p d - 1 ) -p n+1 (-B;p) n 
< A'p n +\-B;p) r 



ln+1- 



Therefore, the inequalities (|3.3|) and (|3.4|) hold for n + 1. Using (|3.4|) we 
reach the following inequality 



^-A^AKAYi-B^U. 



So the limit lim A$ exists. By (|3.2|) . for any 1 < i < d, the lim An' exists. 

n^oo n— >oo 

Thus, we get 



/(*) = /(o) lim 4°, 



which completes the proof. 
Remarks. 
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The condition that f(z) is continuous in \z\ < r can be replaced by the 
assumption that linijv-»oo f{zq N ) exists. 

The above theorem can be easily generalized to multi-variable. 



We now give some examples. The first five examples are adopted from 
the appendix of [27]. 

1. Heine's transformations of 2 4>\ series: 



2Vl 



a, b 



q,z 



{b, az\q) op 
(c, z;g)oo 
(c/b, bz;q) 
(c,z;q) x 
(abz/c;q) c 



2V1 



c/b, z 



az 

abz /c,b 
bz 

c/a, c/b 



;q,c/b 



; q, abz/c 



(3.5) 
(3.6) 
(3.7) 



Let 



We have 



/(*) = 201 



a, b 



-c-q + (qa + qb)z c-qabz 2 
— f\zq) + —, Tt/N )■ 



q(z-l) 



q(z-l) 



(3.8) 



By Theorem 13.11 for \c/q\ < 1, f(z) is uniquely determined by /(0) and the 
recurrence relation (|3.8j) . Let 



(b,az]q)oo 
(c, 2;;g)oo 



c/b, z 



az 



q,b 



Then g(z) satisfies the same recursion as ([3.8)1 . By the g-binomial theorem, 
we have 

c/6 



0(0) = 7 n— i0o 

(c; g)oo 



9,6 



/(o). 



Therefore, ()3.5j) holds for |c/g| < 1. By analytic continuation, ()3.5j) holds for 
all a, b, c, z <G C provided that both sides are convergent. Similar arguments 
can justify (j3.6J) and (|3.7j) . 
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2. Jackson's transformations of 



2<Pl, 292 



and series: 



29l 



a, b 



q,z 



(abz/c; q) 
(bz/c; q) c 
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a, c/b 



q, bz 



-392 



c, az 

a, c/b, 
c, cq/bz 



Q,Q 



(3.9) 
(3.10) 



where (J3.10|) holds provided that the series terminates. 



Let f(z) be the left hand side ()3.9j) . Thus we have the recurrence relation 
(HHI) and lim f{zq N ) = 1. By using the g-Zeilberger algorithm, one can 

N^>oo ' 

verify that the right hand sides of ()3.9|) and ()3.10|) also satisfy the same 
recurrence relation. Moreover, for the summation ()3.10|) the terminating 
condition is required to ensure lim g n ^ = in (j2.1j) . By considering the 

fe— >oo 

limit values, we get the transformation formulas ()3.9j) and ()3.10j) . 

A similar discussion implies the following transformation formula for ter- 
minating 20i series: 



291 



a, b 



q,z 



[c/b, c/a; q) c 



-392 



a, b, abz/c 
abq/c, 



(c/ab, c; q) 

provided that the right hand side summation terminates. 
3. Transformations of 302 series: 



Q,Q 



392 



a, b, c de 

d,e ' ^' abc 
(e/a, de/bc; q) c 

(e, de/abc; q)^ 
{b, de/ab, de/bc; q) c 
(d, e, de/abc; q)^ 



-392 



a,d/b,d/c e 

-q,- 

a 



-392 



d, de/bc 

d/b, e/b, de/abc 



de/ab, de/bc 



q,b 



(3.11) 
(3.12) 
(3.13) 



We take d as the parameter. Let f(d) be the series in (13. lip . We have 
f(0) = l and 



f(d) 



1 + q)ed 2 + {—abc — eb — ea — ec)d + abc + abce/q 
(—ed + abc)(— 1 + d) 

e(— c + dq)(—dq + b)(—dq + a 



+ 



q(-ed + abc)(-l + dq){-l + d) 



f{dq) 

f(dq 2 ). 
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On the other hand, one can verify that both the series in (|3.12j) and (|3.13|) 
have the same limit value and satisfy the same recurrence relation as (|3.11|) . 

4. Sears-Carlitz transformation: 



3^2 



a, b, c aqz 
aq/b, aq/c ' ' be 

(az;q) c 



"5^4 



as, -a2, (ag)2, -(ag) 2 > 
ag/6, ag/c, az, g/2; 



provided that the right hand side terminates. 

Let us take z as the parameter and denote the series by f(z). One can 
verify that both sides have the same limit value lim f(zq N ) = 1 and satisfy 

N—>oo 

the following recurrence relation: 

f(z) = r x {z)f{zq) + r 2 (z)f(zq 2 ) + r 3 (z)f(zq 3 ), 



where 



ab + ac + bc 

r 1 [z) = \-0(z), r 2 {z) 

be 

a 2 

r 3 (z) = — + 0(z). 
be 



-a(b + a + c) 
be 



+ 0(z), 



Note that to comply with the conditions of Theorem 13.11 we only need the 
values ri(0), r 2 (0) and r 3 (0). So we do not give the explicit formulas for 

ri(z), r 2 (z) and r 3 (z). 



5. Transformations of very- well-poised 8 <ft 7 series: 
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a,qah ,—qah ,b,c,d,e, f a 2 q 2 
ah , —a^, aq/b, aq/c, aq/d, aq/e, aq/ f ' ' bcdef 
(aq, aq/ ef, Xq/e, Xq/ /;?)«> 



X 8 7 



(3.14) 



(aq/ e, aq/ f, Xq, Xq/ ef; q)^ 

X,qX^ , —qX^ , Xb/a, Xc/a, Xd/a,e, f _ aq 
A2 , — A2 , aq/b, aq/c, aq/d, Xq/e, Xq/ f ' ' ef 
(aq, b, bcfi/a, bdfi/a, be\i/a, bf\i/a; q)^ 
(aq/c, aq/d, aq/e, aq/ f, \iq, b\i/a; q)^ 

/i, qyh, —q[A, aq/bc, aq/bd, aq/be, aq/bf, bfi/a 
/i^, — /i2", bc\i/a, bdfi/a, be\i/a, bf '/i/a, aq/b ' ' 



(3.15) 



X 8 



(3.16) 



where A = qa /bed and fi = q a /b cdef. 

We choose a, b, f as parameters for the series in ()3.14|) and (|3.15j) and 
denote the series by H(a,b,f). It follows from ()3.12|) that both series have 
the same limit value lim H(aq N ,bq N , fq N ). By computation, one sees that 

they satisfy the following recurrence relation 

H(a, b, f) = r x (a, b, f)H(aq, bq, fq) + r 2 (a, b, f)H(aq 2 , bq 2 , fq 2 ), 

where 

ri (a,b,f) = l + 0(a), r 2 (a,b,f) = 0(a). 

Thus, we have verified the first transformation formula. To prove the second 
transformation formula, we choose a, c, f as the parameters and denote the 
series by H(a,c,f). By computation, the series in ()3.14j) and (|3.1fij) satisfy 
the following recurrence relation 

H(a, c, f) = r x (a, c, f)H(aq, cq, fq) + r 2 (a, c, f)H(aq 2 , cq 2 , fq 2 ). 

Using the transformation formula ()3.13|) . one sees that both sides have the 
same limit value lim H(aq N ,cq N , fq N ). Thus we have obtained the second 

N— »oo 

transformation formula. 



6. A limiting case of Watson's formula. 
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Watson [46] used the following formula to prove the Roger s-Ramanuj an 
identities [30] (see also [27, Section 2.7]): 

y ^^^^r^V ^W^ 72 = (aq; qU V (3.17) 

We choose a as the parameter. Then we can verify that both sides of 
(j3.17|) have the same limit value /(0) = 1 and satisfy the same recurrence 
relation 

/(a) = (1 - aq)f(aq) + aq(l - aq){l - aq 2 )f(aq 2 ). 

Setting a = 1 and a = q in ()3.17j) . we obtain the Rogers- Ramanuj an 
identities by Jacobi's triple product identity: 

OO £.2 oo 

(^^EtV = E(-W® = (g 2 ^ 3 ^ 5 ;g 5 )oo, 

fc=0 W' ^ fc fc=-oo 

and 

OO ^.2_j_^, oo 

Finite forms of the above identities have been proved by Paule [40] by using 
the g-Zeilberger algorithm. 

7. A generalization of Lebesgue's identity. 

The following transformation formula is due to Carlitz [19] (see [2]): 

(x; q) k q® (-a) fc _ (a, x; q)^ (6; q) k x k 
~ (q,bx;q) k {bx-q)^ ^(q,a;q) k 

We choose x as the parameter. Both sides of ()3.18j) have the same limit 
value /(0) = (a; g)^ and satisfy the same recurrence relation: 

/(*) = + /M + ( y + /N 2 )- 

8. Three-Term transformation formulas. 
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Our approach also applies to certain three-term transformation formulas. 
It is sometimes the case that the left hand side of the identity satisfies a 
homogenous recursion, and the two terms on the right hand side satisfy non- 
homogenous recursions respectively but their sum leads to a homogenous 
recurrence relation. 

The first example is 



3^2 



a, b, c 



de 
d,e 

e/b,e/c; q) 



"3<?2 



(e,e/6c; q) 

(d/a, b, c, de/bc; q) 



(d, e, bc/e, de/abc; q) c 



d/a, b, c 
d,bcq/e ' ,q,q 

e/b, e/c, de/abc 
de/bc, eq/bc 



"3^2 



(3.19) 



Let us choose e as the parameter. Then both sides of (j3.19|) have the 
same limit value lim f(eq N ) = 1 and satisfy the same recurrence relation: 

/(e) = r l (e)f(eq) + r 2 (e)f(eq 2 ), 

d 



where 



/x q + d _ . . . 
n(e) = h 0(e), r 2 (e 



Q 



q 



+ 0(e). 



The second example is 



a, qaz, —qa^, b, c, d, e, f 



a 2 q 2 



a * 



2 , —a? ,aq/b, aq/c, aq/d, aq/e, aq/ f ' ' bcdef _ 



-4^3 



aq/bc, d, e, f 
, aq/c, def /a 



(aq, aq/ de, aq/ df, aq/ ef;q) c 
(aq/d, aq/e, aq/ f, aq/def; q)^ 
| (aq, aq/bc, d, e, f, a 2 q 2 /bdef, a 2 q 2 / cdef; q)^ 
(aq/b, aq/c, aq/d, aq/e, aq/ f, a 2 q 2 /bcdef , def /aq; q) c 

aq/de, aq/df, aq/ef, a 2 q 2 /bcdef 
a 2 q 2 /bdef,a 2 q 2 /cdef, aq 2 /def ] q > Q ' 



(3.20) 



We take a,b,f as parameters and denote the series by H(a,b,f). By 
the transformation formula (j3.19|) . we see that both sides of ()3.20|) have the 
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same limit value lim H(aq N ,bq N , fq N ). Moreover, they satisfy the same 

TV— >oo 

recurrence relation: 

H(a, 6, /) = n(a, b, f)H(aq, bq, fq) + r 2 (a, b, f)H(aq 2 , bq 2 , fq 2 ), 

where 

n(a, b,f) = l + 0(a), r 2 (a, b, f) = 0(a). 

9. The Rogers-Fine identity. 

To conclude this paper, we consider a transformation formula that can be 
justified by using nonhomogeneous recurrence relations. This is the Rogers- 
Fine identity [25]: 

( a ; g)fc k _ («, o-zq/b; g)fc(l ~ azq 2k )q k2 - k (bz) k 

h (b,z;q) k (l-zq») " 1 ' > 

We choose z as the parameter. By computation, both sides of ()3.21|) satisfy 
the following recurrence relation: 

(-b + azq)f(zq) + (-zq + q)f(z) = q - b. 

The non-homogenous term g — b appears because g ,o — —Q+b and lim g o k = 

k— >oo ' 

when one implements the g-Zeilberger algorithm. Let d(z) be the difference 
of the two sides of (|3.21|) . Then we have 

d(z) = b - ■ l ' aZqlh d(zq). 
q 1 — z 

Since d(0) = 0, the equation (j3.21|) holds for \z\ < 1 and \b\ < \q\. By analytic 
continuation, it holds for \z\ < 1. 
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